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Abstract. In the framework of the theory of open systems based on completely 
positive quantum dynamical semigroups, we give a description of the continuous- 
variable entanglement for a system consisting of two independent harmonic oscillators 
interacting with a general environment. Using Peres-Simon necessary and sufficient 
criterion for separability of two-mode Gaussian states, we describe the generation and 
evolution of entanglement in terms of the covariance matrix for an arbitrary Gaussian 
input state. For some values of diffusion and dissipation coefficients describing the 
environment, the state keeps for all times its initial type: separable or entangled. In 
^-i other cases, entanglement generation or entanglement collapse (entanglement sudden 

death) take place or even a periodic collapse and revival of entanglement. We show 
that for certain classes of environments the initial state evolves asymptotically to 
■ an entangled equilibrium bipartite state, while for other values of the coefficients 

describing the environment, the asymptotic state is separable. We calculate also the 
logarithmic negativity characterizing the degree of entanglement of the asymptotic 
state. 
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1. Introduction 



X 

When two systems are immersed in an environment, then, besides and at the same 
time with the quantum decoherence phenomenon, the environment can also generate 
a quantum entanglement of the two systems [TJ [2]. In certain circumstances, the 
environment enhances the entanglement and in others it suppresses the entanglement 
and the state describing the two systems becomes separable. The structure of the 
environment may be such that not only the two systems become entangled, but also 
such that the entanglement is maintained for a definite time or a certain amount of 
entanglement survives in the asymptotic long-time regime. 

In the present paper we investigate, in the framework of the theory of open 
systems based on completely positive quantum dynamical semigroups, the dynamics 
of the continuous-variable entanglement for a subsystem composed of two identical 
harmonic oscillators interacting with an environment. We are interested in discussing 
the correlation effect of the environment, therefore we assume that the two systems 
are independent, i.e. they do not interact directly. The initial state of the subsystem 
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is taken of Gaussian form and the evolution under the quantum dynamical semigroup 
assures the preservation in time of the Gaussian form of the state. 

The organizing of the paper is as follows. In Sec. 2 we write and solve the equations 
of motion in the Heisenberg picture for two independent harmonic oscillators interacting 
with a general environment. Then, by using the Peres-Simon necessary and sufficient 
condition for separability of two-mode Gaussian states [31 E], we investigate in Sec. 3 
the dynamics of entanglement for the considered subsystem. In particular, with the help 
of the asymptotic covariance matrix, we determine the behaviour of the entanglement 
in the limit of long times. We show that for certain classes of environments the initial 
state evolves asymptotically to an equilibrium state which is entangled, while for other 
values of the parameters describing the environment, the entanglement is suppressed 
and the asymptotic state is separable. A summary is given in Sec. 4 



2. Equations of motion for two independent harmonic oscillators 

We study the dynamics of the subsystem composed of two identical non-interacting 
oscillators in weak interaction with an environment. In the axiomatic formalism based 
on completely positive quantum dynamical semigroups, the irreversible time evolution 
of an open system is described by the following general quantum Markovian master 
equation for an operator A (Heisenberg representation) [5j[6]: 

= i[H,A(t)] +±.Y,(vH A (t),V j \ + [Vj^m). (1) 

3 

Here, H denotes the Hamiltonian of the open system and the operators Vj, Vj , defined 
on the Hilbert space of H, represent the interaction of the open system with the 
environment. Being interested in the set of Gaussian states, we introduce such quantum 
dynamical semigroups that preserve this set. Therefore H is taken to be a polynomial of 
second degree in the coordinates x, y and momenta p x ,Py of the two quantum oscillators 
and Vj,V] are taken polynomials of first degree in these canonical observables. Then in 
the linear space spanned by the coordinates and momenta there exist only four linearly 
independent operators ^=1,2,3,4 [Z] : 

Vj = a xj p x + a yj p y + b xj x + b yj y, (2) 

where a X j, a y j, b X j, b y j G C. The Hamiltonian H of the two uncoupled identical harmonic 
oscillators of mass m and frequency uo is given by 

H = ^(pl+pl) + ^f( X i + yi). (3) 

The fact that the evolution is given by a dynamical semigroup implies the positivity 
of the following matrix formed by the scalar products of the four vectors a x ,b x ,a y , b y 
whose entries are the components a X j,b X j,a y j,b y j, respectively: 
/ (a x a x ) (a x b x ) (a x a y ) (a x b y ) \ 
(b x a x ) (b x b x ) (b x a y ) (b x b y ) 
(a y a x ) (a y b x ) (a y a y ) (a y b y ) 
\(b y a x ) (b y b x ) (b y a y ) (b y b y )J 



2 



(4) 
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We take this matrix of the following form, where all coefficients D xx , D xpx} ... and A are 
real quantities (we put, for simplicity, h = 1): 

/ D xx -D xPx -i\/2 D 

-£>™„ + iA/2 



D 



'xp, 

D 



-D 



X Py 



j-y 



PxPx 

Dyp x 



-D 
D 



yp.r 



D 



PxPy 



illl 



-D 



V 



D X py Dp x p y 



D yp + iA/2 



ypy 
D 



iA/2 



(5) 



Py'Py 



J 



It follows that the principal minors of this matrix are positive or zero. From the Cauchy- 
Schwarz inequality the following relations for the coefficients defined in Eq. ([5]) hold: 



D XX D PxPx 

D xx Dyy 
DyyD px p x 



D 



■IP:, 



A 2 

— DyyD py p y 



X 2 

D 2 >—, 

ypy — 4 ' 



DL > o, /;, ,./; 



xy 



PyPy 



Dl 



D 



HI'., 



> 0, Dp xPx Dp yPy 



■ry > °> 

- D lxPy > 0- 



(6) 



(7) 



The matrix of the coefficients ([5]) can be conveniently written as 

( Ci c 3 \ 

in terms of 2 x 2 matrices C\ = C\\ C 2 = C 2 and C3. This decomposition has a direct 
physical interpretation: the elements containing the diagonal contributions C\ and C 2 
represent diffusion and dissipation coefficients corresponding to the first, respectively the 
second, system in absence of the other, while the elements in C3 represent environment 
generated couplings between the two, initially independent, oscillators. 
We introduce the following 4x4 bimodal covariance matrix: 



a(t) 



( 


0~xx 


Oxpx 


0~xy 


0~XPy 


\ 




a xp x 


<J P X Px 


a ypx 


OpxPy 






0~xy 


°ypx 


a yy 


a yp y 




\ 


°xpy 


OpxPy 


a yp y 


a PyPy 


) 



(8) 



with the correlations of operators A\ and A 2 , defined by using the density operator p 
describing the initial state of the quantum system, as follows: 

1 



i Tr(p(A 1 A 2 + A 2 A l ){t)) - Ti(pA 1 (t))Ti(pA 2 (t)). 

By direct calculation we obtain [7] (T denotes the transposed matrix): 
da it) 



dt 



Ya(t) + a(t)Y l + 2D, 



where 



Y 



D 



( ~ x 


1/m 





\ 


—muj 2 


-A 














-A 


1/m 


V 





—mu 2 


-A ) 


^ -D X x 




D X y 


D X p y ^ 


Dxp x 


Dp x p x 


Dyp* 


Dp xP y 


D X y 


Dypx 


Dyy 


Dypy 


\ D X py 


DpxPy 


Dypy 


DpyPy / 



(9) 



(10) 



(12) 
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= ( d). ( 15 ) 



The time-dependent solution of Eq. (fill is given by [7] 

a(t) = M(t)O(0) - a(oo))M T (t) + a(oo), (13) 

where the matrix M(t) = exp(Yt) has to fulfill the condition lim^oo M(t) = 0. In order 
that this limit exists, Y must only have eigenvalues with negative real parts. The values 
at infinity are obtained from the equation [7J 

YV(oo) + a(oo)Y T = -2D. (14) 
3. Dynamics of entanglement 

The two-mode Gaussian state is entirely specified by its covariance matrix (JED, which is 
a real, symmetric and positive matrix with the following block structure: 

A C 

, ° T B , 

where A, B and C are 2x2 matrices. Their entries are correlations of the canonical 
operators x,y,p x ,p y , A and B denote the symmetric covariance matrices for the 
individual reduced one-mode states, while the matrix C contains the cross-correlations 
between modes. The elements of the covariance matrix depend on Y and D and 
can be calculated from Eqs. f|T3|) . ffT4"j) . Since the two oscillators are identical, it is 
natural to consider environments for which the two diagonal submatrices in Eq. (J7j) 
are equal, C\ = C2, and the matrix C3 is symmetric, so that in the following we take 
D xx = D yy ,D xPx = D yPy ,D PxPx = D PyPy ,D xPy = D yPx . Then both unimodal covariance 
matrices are equal, A = B, and the entanglement matrix C is symmetric. 

3.1. Asymptotic entanglement 

First we analyze the existence of entanglement in the limit of large times. With the 
chosen coefficients, we obtain from Eq. ( fl4l) the following elements of the asymptotic 
entanglement matrix C(oo): 

/ v _ ™ 2 (2\ 2 + u 2 )D xy + 2m\D xpy + D PxPy 
a ^°°>- 2m 2 X(X 2 + u 2 ) ' 1 j 

( , ( s -m 2 u 2 D xy + 2mXD xp +D PxP 
<W«0 = ^(00) = 2m(A2 + L , 2) » ( 17 ) 

f s m 2 u 4 D xy - 2mu 2 XD xpy + (2A 2 + uj 2 )D PxPy 
^ (00) = 2X(X 2 + u 2 ) • (18) 

The elements of matrices A(oo) and B(oo) are obtained by putting x = y in the previous 
expressions. We calculate the determinant of the entanglement matrix and obtain: 

1 



det C(oo) 



4X 2 (X 2 + lu 2 ) 



x \{mu 2 D xy + ±-D PxPy ) 2 + AX 2 (D xy D PxPy - D 2 )]. (19) 



in 
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It is interesting that the general theory of open quantum systems allows couplings 
via the environment between uncoupled oscillators. According to the definitions of 
the environment parameters, the diffusion coefficients above can be different from zero 
and therefore can simulate an interaction between the uncoupled oscillators. Indeed, 
Gaussian states with det C > are separable states, but for det C < it may be possible 
that the asymptotic equilibrium states are entangled, as will be shown next. 

On general grounds, one expects that the effects of decoherence, counteracting 
entanglement production, be dominant in the long-time regime, so that no quantum 
correlation (entanglement) is expected to be left at infinity. Nevertheless, there 
are situations in which the environment allows the presence of entangled asymptotic 
equilibrium states. In order to investigate whether an external environment can actually 
entangle the two independent systems, we use the partial transposition criterion [3j|4]: 
a state is entangled if and only if the operation of partial transposition does not preserve 
its positivity. For the particular case of Gaussian states, Simon [I] obtained the following 
necessary and sufficient criterion for separability: S > 0, where 

S = detA det B 

+(- - | det C\f - Tr[AJCJBJC T J] - -(det A + det B) (20) 
and J is the 2x2 symplectic matrix 

-(-,:)■ 

In the following we consider an environment characterized by diffusion coefficients of 
the form m 2 uj 2 D xx = D PxPx , D xPx = 0, m 2 uo 2 D xy = D PxPy . This corresponds to the case 
when the asymptotic state is a Gibbs state [6]. Then Simon expression (1201) takes the 
following form in the limit of large times: 



X 2 X 2 + oo 2 A) 

a ^ D XX D xpy / nn\ 

\ 2 (\ 2 + UJ 2 ) ' 1 ; 

For environments characterized by such coefficients that the expression S(oo) (|22p is 
strictly negative, the asymptotic final state is entangled. In particular, for D xy = we 
obtain that S(oo) < 0, i.e. the asymptotic final state is entangled, for the following 
range of values of the coefficient D xpy characterizing the environment jH [9] : 

muD xx 1 D xpy muD xx 1 

2 < 7F+^ < ^^ + 2' {26) 

where the diffusion coefficient D xx satisfies the condition muD xx /X > 1/2, equivalent 
with the unimodal uncertainty relation. We remind that, according to inequalities ([6]), 
the coefficients have to fulfill also the constraint D xx > D xpy . If the coefficients do not 
fulfil the inequalities (1231) . then S(oo) > and the asymptotic state of the considered 
system is separable. 
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Figure 1. Asymptotic logarithmic negativity E(oo) versus environment coefficients D 
and d. We have taken m = u> = Ti = 1 and A = 0.2. According to inequalities ([6]), only 
the part of the plot corresponding to the constraint D > d has a physical meaning. 



3.2. Asymptotic logarithmic negativity 



We apply the measure of entanglement based on negative eigenvalues of the partial 
transpose of the subsystem density matrix. For a Gaussian density operator, the 
negativity is completely defined by the symplectic spectrum of the partial transpose 
of the covariance matrix. The logarithmic negativity E = — | log 2 [4/(cr)] determines 
the strength of entanglement for E > 0. If E < 0, then the state is separable. Here 



f(a) = -(det A + det B) - det C 



; (det A + detB) -detC 



1/2 



det a 



In our case the asymptotic logarithmic negativity has the form 

mujD TT D, r 



Moo 



log 2 



'X P y 



(24) 



(25) 



It depends only on the diffusion and dissipation coefficients characterizing the 
environment and does not depend on the initial Gaussian state. As an example, in 
Figure 1 we represent the asymptotic logarithmic negativity E(po) versus diffusion 
coefficients D xx = D and D xp = d. We notice that for some range of these coefficients, 
E{oo) takes strictly positive values, measuring the degree of entanglement of the 
corresponding asymptotic entangled states. Again, the coefficients have to fulfill the 
constraint D xx > D xPy , so that in Figure 1 only the part of the plot corresponding to 
D > d has a physical meaning for our model. 
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-0 . 5 

Figure 2. Simon separability function S versus time t and environment coefficient d, 
for A = 0.5, D = 0.4 and a separable initial Gaussian state with initial correlations 
cr xx (Q) = 1, ct PxPx (0) = 1/2, crxp x (0) = cr xy {0) = <r Px p y (fy = &xp y (0) = 0. We have 
taken m = u> = h = 1. 

5.5. Time evolution of entanglement 

In order to describe the dynamics of entanglement, we have to analyze the time evolution 
of the Simon function ( 1201) . We consider separately the two cases, according to the type 
of the initial Gaussian state: separable or entangled. 

1) To illustrate a possible generation of the entanglement, we represent in Figure 
2 the function S(t) versus time and diffusion coefficient D xp = d for a separable initial 
Gaussian state with initial correlations a xx (0) = 1, o~ PxPx (0) = 1/2, cr xpx (0) = o~ xy (0) = 

= <7 xp y (0) = 0. We notice that, according to Peres-Simon criterion, for relatively 
small absolute values of the coefficient d, the initial separable state remains separable 
for all times, while for larger absolute values of d, at some finite moment of time the 
state becomes entangled. In some cases the entanglement is only temporarily generated 
and the state becomes again separable after a certain moment of time. In other cases, 
namely for even larger absolute values of the coefficient d, the generated entangled state 
remains entangled forever, including the asymptotic final state. 

2) The evolution of an entangled initial state is illustrated in Figure 3, where 
we represent the function S(t) versus time and diffusion coefficient d for an initial 
entangled Gaussian state with initial correlations a xx (0) = 1, cr PxPx (0) = 1/2, o~ xpx (0) = 
0, o~ xy (0) = 1/2, cr PxPy (0) = —1/2, a xp (0) = 0. We notice that for relatively large 
absolute values of the coefficient d, the initial entangled state remains entangled for all 
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-0 . 5 

Figure 3. Same as in Figure 2, for an entangled initial Gaussian state with initial 
correlations a xx (0) = 1, cr PxPx (0) = 1/2, a xPx (0) = 0, a xy (0) = 1/2, o PxPv {0) = 
-1/2, a xPy (0) = 0. 

times, while for smaller absolute values of d, at some finite moment of time the state 
becomes separable. This is the well-known phenomenon of entanglement sudden death. 
Depending on the values of the coefficient d, it is also possible to have a repeated collapse 
and revival of the entanglement. 

4. Summary 

In the framework of the theory of open quantum systems based on completely 
positive quantum dynamical semigroups, we investigated the existence of the quantum 
entanglement for a subsystem composed of two uncoupled identical harmonic oscillators 
interacting with a general environment. By using the Peres-Simon necessary and 
sufficient condition for separability of two-mode Gaussian states, we have described 
the generation and evolution of entanglement in terms of the covariance matrix for 
an arbitrary Gaussian input state. For some values of diffusion and dissipation 
coefficients describing the environment, the state keeps for all times its initial type: 
separable or entangled. In other cases, entanglement generation or entanglement 
collapse (entanglement sudden death) take place or even one can notice a repeated 
collapse and revival of entanglement. We have also shown that, independent of the 
type of the initial state, for certain classes of environments the initial state evolves 
asymptotically to an equilibrium state which is entangled, while for other values of 
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the coefficients describing the environment, the asymptotic state is separable. We 
determined also the logarithmic negativity characterizing the degree of entanglement 
of the asymptotic state. The existence of quantum correlations between the two 
considered harmonic oscillators interacting with a common environment is the result 
of the competition between entanglement and quantum decoherence. 

Due to the increased interest manifested towards the continuous-variable approach 
to quantum information theory, these results, in particular the possibility of maintaining 
a bipartite entanglement in a diffusive-dissipative environment for asymptotic long 
times, might be useful in controlling the entanglement in open systems and also for 
applications in quantum information processing and communication. 
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